PROJECTION METHODS FOR SOME CONSTRAINED SYSTEMS 
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Abstract. This article is concerned with a geometric tool given by a pair of projector 
operators defined by almost product structures on finite dimensional manifolds, polarized 
by a distribution of constant rank and also endowed with some geometric structures 
O , (Riemann, resp. Poisson, resp. symplectic). The work is motivated by non-holonomic 

' and sub-Riemannian geometry of mechanical systems on finite dimensional manifolds. 

, Two examples are given. 



o 



1. Introduction 



, In general, a classical constrained mechanical system consists in three basic ingredients: 

an n-dimensional configuration manifold W, a polarization TD on W, which is to say a 
distribution P : z G W — > C T^W and an auxiliary geometric structure. V is hereafter 
supposed smooth, of constant rank (dim = constant, for all z), and the Lie algebra of 
^ ! vector fields taken as sections of T) over W span the tangent space of W at each point 

ly^ \ (thus V is non-integrable, i.e. non-involutive in the Frobenius sense). We recall that a 

00 ■ curve (p(t) in W satisfies the constraints (is polarized) if (p(t) G ^^(p(t)- 

With these assumptions, we may outline some interesting directions on polarized sys- 
tems: suppose first that the auxiliary geometric structure is given by a Riemannian metric 
tensor g such that TW splits as TW = D©I>-'-, where T)-^ is the g-orthogonal complement 
of V in TW: g(v, w) 0,v 6 2^z)W G P^, for all z G W (we use the same notation for 
the distribution and the corresponding sub-bundle of TW). Sometimes, the polarization 
^ ! D is said horizontal, and Vaisman [IS], [IH], following Reinhardt |1D], coined the triple 

\ (W, g) a Riemannian almost foliated manifold. If V is Frobenius integrable, (W, "D, g) 

is said a Riemannian foliated manifold. 

In classical mechanics, the triple (W, P, g) defines two different important mathematical 
structures: 

(1) Non-holonomic (NH) mechanics by assuming that the trajectories satisfy D'Alem- 
bert's principle of virtual work: the constraining force must be perpendicular to 
the horizontal subspace, since it does not produce work (see §4.1j see also ref. [2Tj . 
p. 85 or ref. ^34]). 

(2) Vakonomic (YAK) mechanics by assuming that the trajectories do not obey D Alem- 
bert principle and satisfy a Lagrange variational principle (see Arnol'd et. al. [2] 
for further information). 

Unless the distribution V is Frobenius integrable, VAK mechanics gives different geodesic 
equations from NH mechanics, and the comparison between these structures was elluci- 
dated in ref. [TB] (previous works on the subject are, for instance, references P2] and [T7]). 
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Recall from ref. |18] that YAK mechanics is related to the so called sub-Riemannian (SR) 
or Carnot-Caratheodory geometry (see the book edited by A. Bellaiche and J-J. Risler 
[5], R. Montgomery [33], I. Kupka [25] for further details in SR geometry, and Koiller et 
al. [23] for an example in sub-Riemannian Lagrange mechanics). 

Now, the splitting TW — V (B T)-^ means that W admits an almost product structure 
a.p.s., for brevity, i.e, W is endowed with a tensor field T of type (1,1), taken as a vector 
valued one form on the tangent bundle P: TW TW, of involutive character, i.e, such 
that = id (Schouten [51], Nickerson-Spencer p5] and Walker [51]). In fact, there are 
defined the following bundle maps: a P^- valued one form P: TW V-^, (a tensor field 
of type (1,1) on W), such that P is a projection operator onto P"*": P o P = = P, and 
a complementary Q = id — P: TW T). Then the pair (P, Q) defines P = Q — P, with 
eigenvalues 1 and —1, and associated eigenspaces V and P-*-, respectively. We remark 
that the interplay of a.p.s. with covariant derivatives is one of the themes of Hermann's 
book [21] (see also ^). 

On other hand, let us suppose that the polarized manifold (W, P, g) is endowed with 
an a.p.s. F, and in addition the eigenvector bundle corresponding to the eigenvalue 1 is 
precisely P^, at each point z of W. Then we have the projectors 

Q = (i/2)(id + r),p = (i/2)(id-r), 

and T^W = V^® P^, Vz G W, where V = imP, image of P in TW. Therefore, an 
alternative case may be considered, where the subspaces C T^W are hereafter supposed 
of complementary constant dimension, everywhere transversal to those of V, for all z G W. 
For simplicity, we shall call this decomposition by oblique. 

One more polarized situation is illustrated by Poisson manifolds. Consider the pair 
( W, TT) , where TT is a twice contravariant skew-symmetric tensor field, verifying [TT, TT] = 
0, where [, ] are the Schouten brackets (see ref. [32], [H] for further details). The tensor TT 
is called Poisson's tensor field or Poisson structure, and ( W, TT) a Poisson manifold. The 
Poisson structure induces a bundle morphism jj : T*W TW such that |3 jj (oc) = n(a, (3), 
where a and (3 are one forms on W. In particular, dg jj (df) = TT(df, dg) = {f, g} is the 
well known Poisson bracket of f , g G C°°(W), the space of C°° functions on W. 

Let us suppose that the Poisson bivectors are of constant rank < n. Then the char- 
acteristic distribution «S: z G W i— > tjz^(T*W) is differentiable and completely integrable, 
and defines a foliation of W such that each leaf S is endowed with a unique symplectic 
structure and the tangent space T^S through each point z is |jz(T*W). Therefore, we 
may also consider the case where the a.p.s. is defined by a polarization (W, 5, TT), such 
that T^W = \S © S^, with the second factor being a complementary distribution of con- 
stant dimension. In fact, these comments motivate a study of the inter-relation of almost 
product geometry with some fields of classical mechanics. 

Indeed, the purpose of this article is to retake this subject for the cases where the 
projectors can be modeled by a Riemannian or a Poisson (resp. symplectic) structure 
on a finite dimensional manifold. This means that we will be in the context of the 
determination of appropriate projectors defined by these geometric objects. We propose 
as a first task to re-examine the relationship which exists between the almost product 
geometry with a (non-degenerate) Riemannian structure. We adopt the Pfaffian view 
point of ref. [33], and we assume that (W, g), endowed with a set of n — m-linearly 
independent one-forms w*, such that — defines a non-involutive distribution of 
constant dimension, with a complementary integrable distribution, i.e, W is a foliated 
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Riemannian manifold. We shall express the metric in terms of a local coframe {dz'', cu"} 
on the configuration manifold. Therefore W admits an "oblique a.p.s.". Naturally we 
may orthogonalize the a.p.s., i.e, the a.p.s. is defined by a cobasis of orthogonal covectors 
with respect to the given metric. However, as we shall see in the examples, we may work 
directly with the oblique situation to obtain the projected dynamical equations, avoiding 
the application of the Gram-Schmidt procedure. This will be not only convenient for 
matrix calculations, but also to set up other studies like the equivalence problem (ref. 
[53!, Particular cases are the sub- Riemannian geodesic problem and the Vakonomic 

variational approach. 

Next, as a second task, we consider an involutive polarization on a Poisson manifold, 
defined by its integrable symplectic foliation. We apply the projector method to obtain 
the Poisson structure for the transverse manifold of the symplectic leaves. This could be 
seen as an application of the previous study in the sense that the metric is replaced by 
a Poisson structure, and the integrable distribution is the characteristic distribution. We 
profit the occasion to give the Dirac formula for constraint manifolds in the transverse 
situation (see ref. [32] for the symplectic submanifold case; the technique used is the usual, 
but we have searched the literature and have not found it for the transverse situation). 
The reader will find more about Dirac mechanics in non-holonomic contexts in references 
[B] and [S]. 

Finally, we would like to stress here that the projector method is an adequate tool to 
treat some variational problems in which the extremal curves as well as the comparison 
curves (associated with a given Lagrangian function) are required to fulfill conditional 
equations (the constraints). As it is well known, the oldest problem of this type was solved 
by Pappus in the third century A.D. These variational problems are called The Problem 
of Lagrange, who first formulated the problem clearly [11]. The projector method is 
based on the (orthogonal) decomposition of the virtual displacement, which make it right 
for dealing with a generalized form of D'Alembert's Principle from which the equation 
of motion of mechanical systems are usually derived (see ref. [IB] and also the book of 
Arnol'd, ref. [1], p. 91 —95). We remark that, when the method is applied to the Problem 
of Lagrange, neither Lagrange's multipliers nor elimination of coordinates is required to 
obtain the equations of motion. In this way, this approach, is appropriate to implement 
the canonical quantization of constrained systems, because the ambiguities introduced by 
the Lagrange's multipliers are eliminated. It is also appropriate for setting up computer 
calculations for large multibody systems which appears in control problem of mechanical 
systems and robotics (for example, see ref. [13] for a computer use of IEEE Scheme 
Programming in Mechanics) . 

This paper is structured in three sections. In section 2 we examine some intrinsic and 
local properties of the projector method on a Riemannian manifold. We begin with some 
intrinsic considerations using the so-called musical bundle morphisms jjg, bg, induced by 
the Riemannian metric g. Next we suppose that D is integrable, characterized by a set of 
k-linearly independent one-forms. We express the pair of bundle projectors in terms of this 
set, we obtain the local expressions for the corresponding bundle projection morphisms, 
and then we give a local description in terms of the Riemannian metric. 

In section 3, we replace the geometric structure g by FT, and we study the role of the 
bundle projectors to obtain the Poisson structure for the transverse manifold of the sym- 
plectic leaves. We shall return to this situation in the second example of the last section 4. 
Indeed, this section is only devoted to applications of the method. We start subsection 4.1 
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with a brief review on D'Alembert's Principle and then we study the so-called Chaplygin- 
Caratheodory sleigh, a prototype of a non-holonomic constrained system. Next, in sub- 
section 4.2, we examine a system consisting of a free particle in IR^, subjected to the 
non-holonomic contact form w = dz — ydx. From the projector's viewpoint one obtains 
a very well known non-holonomic Lie algebra, the Heisenberg algebra, a fundamental ex- 
ample in sub-Riemannian geometry. We conclude the example with a study of the motion 
of the particle in the phase space using the underlying Poisson structure. 

The following convention will be adopted: capital roman letters I, J, K, etc. run from 1 
to n. Lower case roman characters a, b, c run from 1 to m, representing the constraint 
distribution. Greek characters a, |3, y, etc., run from 1 to n— m. Summation over repeated 
indices is assumed unless otherwise stated. By a differentiable manifold, we shall mean 
C°°, connected, separable and Hausdorff. 

2. Polarized Riemannian manifolds 

2.1. An intrinsic relation. Let (W, P, g) be a n-dimensional Riemannian manifold, 
supposed endowed with an a.p.s. F, compatible with V in the sense that the eigenvector 
bundle corresponding to the eigenvalue — 1 is precisely X'z, at each point z of W. Let us 
take the corresponding bundle projections P, Q = Id — P so that Q : TW — > V projects 
onto V, and P: TW projects onto a complementary distribution V^. 

Now, the tensor g defines a bundle isomorphism jjg : T*W TW with inverse denoted 
by bg: TW T*W (the so-called musical morphisms), defined respectively by 

m) = g.(*, •] = Z4,, bg(z) = g*(., z) 4>z. 

As g is symmetric, one has [Ig = ttg, where T*W — > TW is the adjoint operator 
= (1) o ttg. So, if we set ftq = Q o ftg, ftp = P o ftg, then 

ttg = Qottg = (QTottg = ttgoQ* = fQ 
ftp = Poftg = (PTottg = ttgoP* = ft*, 

and it follows that 

(1) gJQ^(4)),il.) = •ilj(ttgoQ*)4) = il.(Qoftg)4) = Q'^(i|;)ttg4) = g.(4),Q*(i|;)] 

g*(Q'(*),-^] = l^(ttgoQ*)4)=lIj(ttgo(Q^)2)4)=v|;(QottgoQ*)cl) 

(2) = Qn^P]ttgQ*(*]=g.(Q*(*],Qn^l^]). 
Also, 

g,(P*((|)),P*(il))) = g.(l(4))-Q*(4)),l(i|;)-Q^(ct))) 

= g.((t),ip)-2g.(Qn4)),tP)+g.(Q*((l)),Q*(iP)) 

(3) = g*(*,^^)-g*(Q*(*),Q*(^l>)). 

If we set 

gQ = g*(Q*(«), Q*(»)), gp = g*(P*(«), P*(»)) 

then g^ = gp+gQ, as it would be expected, i.e. the almost product structure T* = P*— Q* 
is such that g*(r*(4)), r*(i|;)) = g*(4>)''J^) for aU ({), ij>. Obviously, 

(4) ftg=ftQ+ttp. 
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Remark 1. We observe that if we replace the tensor g by a Poisson tensor FT or by a 
symplectic structure H then we have a similar result if the above compatibility assumption 
on r is assumed (see also p. [8]). □ 

2.2. Local expressions for the case of a foliated manifold. Let us suppose now 
that V is a completely integrable distribution on W, of constant rank n — m, and so 
(W, g, r) is a Riemannian foliated almost product manifold. The distribution defines 
a foliation on W, denoted also by the same symbol V, to simplify things. The tangent 
bundle TV is the vector sub-bundle of TW such that "D^ = Tz£z for any leaf £ of V and 
any z ^ £. 

Let us consider the splitting T^W = © P^, z G W, where is a distribution of 
subspaces C Jz^, of complementary constant dimension, everywhere transversal to 
those of T>. To look for the local expressions, let U be a neighborhood of z G W, so 
that (due to the integrability of V) the leaf is locally given by equations z'^ = 0. Thus 
we write = span{Yc^ = 9/9z*}. Let be a set of linearly independent vectors with 
= span{Ya}. Furthermore, writing Yq = (9/9z^) in the coordinate basis one obtains 
9/9z" — 9/9z*, as a new basis, for suitable functions F*(z) on W. Then we define the 

v.: 



projectors P: T^W - 
(5) P(Z] 

with z = z^fg/gz'^) 

(6) 



= z^(^-i 

f Z«(9/9z*). 



P: T,W 

, 9 , 



Q(Z) 



9z« 

In matrix notation: 



(Z^ + F^Z''] 



9z'^ 



id 
-F* 



To simplify the notation we set heretofore Xa 
illustrates the situation: 




F« id 

d/dz"- — F*9/9z*. The following figure 




(Q-P) (Zz 



Let = dz* + F"(z) dz'^ be the set of independent 1 -forms such that {dz'',w"} is the 
corresponding cobasis for the cotangent space T*W. Then 

T*W = [VIT ® i'DzY = span{dz"} © spanlw*"}, 
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and we may set 

P = dz'^ ® X„ Q = w« ® Y„. 

We observe that locally (see (jUand ([5])), 

Next we shall study these projectors in terms of the Riemann matric g — gijdz^ (8) dzK 
To do this we recall that g has the following expressions (see ref. |45j): in the basis 
{dz'', w«} 

(7) g^ = (g.b - 2 r« + g„p T^jdz^^ ® dz'' 

+2 (ga« - gap r^^) dz'^ ® + g„p w« ® 

= Fab dz" ® dz^ + Fa« dz" ® w« + Gocp w« ® 
and in the basis {X^, Y^}, 

g, = g"^Xa®Xb+2(g-« + g-^r«)Xa®Y„ 

^(gap ^ 2g- rp + g'^^ r« rp) y, ® Yp 

•ilf G'^' Xa ® Xb + G"* Xa ® Y„ + G«P Y„ Yp. 

Let = g^rfw", •). Then w'^(£,*) = g*(w", w'^) = G*'^. As the matrix with entries 
gap jg invertible, let G^p be the set of functions which are the entries of the inverse. We 
define the tensor 

(8) q = g*(.,e)GapgJw^.). 

Then one obtains 

(9) q = G„pW*®£,P, 
and so q(£,P) - (G^" G^y) U- As 

= [g^(^nG„pg.(w^.)][g^(.,e)G^eg.(w^.)] 
= g^(., Gap [g.(wP, .) g^., e)] G^e g.(w^ .) 
= g'^(^e)GapGP^G^eg.(w^.) = q, 

q is a projector onto the space spanned by the £,°"s, with complementary projector p. 
Let us suppose now that gaa— gap vanishes in ([7]). Then and VI are g-orthogonal, 
= X>-^, g*(Xa,Ya) = and so g* admits the following diagonal form with respect to 

{dz^w-}, 

g^ = Gab dz" ® dz^ + Gap ® wP, Gab = g*(Xa, Xb), G^p = g^^lYa, Yp) 

(or gj, = G'^'' Xq (8) Xb + G"'^ Ya (8) Yp, where G"'', resp., G*'^ are the entries of the inverse 
matrix of (Gab, resp., (Gap)). As = g^(w", •) = G""*' Yy, then it is easily verified that 

(10) q = w« ® Ya = Q 

and so Im q = span{Ya} (and obviously p P), i.e, these projectors have the same local 
matrix expression given by ([6]). 
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3. Projectors and Transverse Poisson structures 

In this section we first consider a polarization on a Poisson manifold, defined by its 
integrable symplectic foliation, and the local expression of the corresponding projectors 
in terms of the symplectic form. In the second part we study the role of the bundle 
projector TW|M TM, where M C W is a given transversal (holonomic) manifold of 
codimension = dimension of the symplectic leaf (see the Introduction), in the process of 
reduction of a Poisson manifold (W, TTw)- We use the projector method to obtain the 
local expression of the corresponding induced Poisson structure (MjTTm)- 

Throughout this section we follow the following convention: local coordinates are now 
denoted by z = (z'^, z^), the characters a, b,c running from 1 to m and u, v from 1 to 
n — m = k. Greek characters a, |3,y, etc., are used for differential forms on manifolds. 

3.1. Projectors for the symplectic foliation. We assume that the smooth manifold 
W is endowed with a Poisson structure - hereafter denoted by TTw - of constant rank 
m < n, which induces a bundle morphism jj: T*W TW such that |3 {t (a) = TTw((^> P), 
where a and (3 are one forms on W. 

Let S be the unique symplectic leaf of the characteristic distribution z i-^ ttz(T*W) 
going through z G W. Thus S is obviously a Poisson submanifold. Now, as the rank of 
TTw is constant, the symplectic leaves are of constant dimension = m, and we may choose 
a decomposition 

\w = T,s e s^, 

with T^S = span{Ya}, where Ya = d/dz"", and = span{Xu = 9/9z^ - X.l=^ K Ya), for a 
local coordinate chart (U, (z'', z^)). The Poisson bivector is locally TTw = j "^^^ A ^h- 

So, from this local expression one obtains an almost symplectic form & (see Vaisman 
[H], p. 37) expressed as O = ^AabW^ Aw'', with w'' being the dual form of Yq, and 
tt^'^Aqc = Sc (recall that n"-^ = Therefore one obtains the following projectors : 

the first, taking into account that 

0(., Y-) = A,b w^(YJw'' = A,b S^*^ = A,b w^ 

is given by (compare with Q) 

e(.,Y-)Aab^w(w^.) 

(11) 

and the second (compare with ffTOl) ). 

e(.,£,'')A,b^w(w^•] = w'^ Aab^,'' =Aab7r^'w'^®Yc 
(12) = w^®Ya = G(•,r)Aab^w(dz^.). 

3.2. The transverse holonomic case. Let us suppose now that M is a submanifold of 
the polarized Poisson manifold (W, iS, TTw), dJid denote by 

Ann T^M = {cc^ e T^W; aJT^M) = 0}, 

the annihilator of TM C TW|jvi in T*W. Suppose that 

(a) T^M n Hz Ann T^M = {0}, (b) ker jj^ n Ann T^M = {0}, Vz G M 



= e(., Y-) A,b e = e(., Y")A,b 7r^^ Yc 
= A,bW^®Ya = e(«,Y'^)A,bnw(dz\.), 
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which are equivalent to the condition T^W = TjlVl © jj^ Ann T^M. (see ref. |T9j, p. 126, for a 
proof). Then Weinstein ([52], P- 529) proved that M. is endowed with a Poisson structure 
TTm, defined by the composition 

(13) riU ^ rW|M ^ TW|M ^ TM 

where ptm is the bundle projection along HzfAnnT^M) onto TM and p\f^ is its adjoint 
(Proposition 1.4 of ref . [52], p. 529/530). Now, it can be shown that (see ref. [H], p. 39) 

- the assumption (a) is equivalent to the statement that T^M fl T^S is a symplectic 
subspace of T^S, and 

- the assumption (b) is equivalent to the statement that M is transversal to the 
symplectic leaf S passing through z. 

Therefore \W = \M + \S, and: 

- the Poisson tensor of W is the product of the Poisson tensor TTs by Hm, the 
so-called transverse Poisson structure to S at z (TTs is induced by the symplectic 
structure of S), 

- the transverse Poisson structure of M. may be computed via the Dirac's bracket 
formula if further assumptions are made (see refs. [31], [32], Prop. 8.5.1, p. 226, 
[36j . Prop. 2, p. 88, [15] or even ref. [38] )• In such a case, one has the so-called 
Dirac's theory of second class constraints. 

Let us compute the transverse Poisson structure from our viewpoint, but supposing 
that the dimension of M, going through z, is the codimension in W of the corresponding 
symplectic leaf S of the foliation and so 

\W = TM ® T^S, 

(if M is transversal to S so that T^M fl jlz Ann T^M ^ {0} is a distribution of constant rank 
then we shall need more assumptions - see Vaisman [44j, for instance). 
Consider the following composition map 

tls = ttop|s:rS^TW|s. 

If = I — pyg is the complementary projector and if we set 

ttM = ttop^^:rM-^TW|s, 

then obviously tts + jJivi = tl and for all forms a, (3 one has 

atls(3 + (xUmP = (xtJP = nw(a, |3). 

This expression suggests that (x|ts|3 (resp. ajjivip) is a good candidate for the Poisson 
tensor Hs (resp. TTm)- In fact, if we suppose that H o p^g = pjg o jj then it is easy to 
show that ws (tts(c>c), tts(|3)) = cxttsP, where wg is the symplectic structure on S, and so, 
from Proposition 3.2 of ref. [38], one has TTwfPTs i'^^i Pts (P)) = c<.tlsP- 

Let us set TTs = Pts (Hw), that is, 

ns(a, |3) =nw(p^s'^)PTsP) = a^sp. 

Then 

Hm = Ptm (Hw) = (I - Pts) (TTw) = TTw - Tls 

gives the complementary relation TTmIoc, (3) = (xjjjvip. 
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The above tensor TTm is a Poisson tensor, as a consequence of the symplectic structure 
of S. So, let us see this in terms of local conditions related to the symplectic manifold S: 
let Zo G U, an open subset of W and (x\ • • • ,x^^) local coordinates for S, used to define 

MnU ={zG U c W-.x'^lz) =0,a = 1,--- ,2s}, 

transversal to S, through zo . Then, 

(14) ws(X%X^) =ns(dx%dx^) ={x%xV = X^(x'^) = dx^iX^) = A''\ 

Here X'^ is the Hamiltonian vector field associated to dx'', by the bundle homomorphism 
tts defined by wg. For what follows, we shall denote by A^b the entries of the inverse 
matrix of (A''^). 

Let f e C°°(M) be such that {x^f}w = X"(f) = df(X^) = and f e C°°(W) an 
extension of f to a neighborhood in W, written as 

(15) f = f + UbX^ 

so that the action of X" on both sides of ( ITSl) gives Ub = — X"(f)Aab, and so 

f = f + X'^if) AabX^ = (I - X^ Aba X'^]? . 

Now, the tensor field t = (I — x^'AbaX'') is a projector. To see this it is sufficient 
to show that a = x'^AbaX'' is the complementary projector. Indeed, ct(x'^) = x'^ and as 
(T(g) =X'^(g)AabX^ then 

0^(9) = X^(g)AabCT(x") 

= X^(g] [Aab A^^ Aed x"] = X-^lg) A.b 6"d 

= X(g)AabX^ = cT(g), 

and obviously t(x") =0, Tocr = croT = 0. Thus the projector pjg = T*S T*W|jvi is 
taken as Pjg = dx'^ Aab X^ (or Aab dx'' X^) and so from (IT^ we have, 

p^s (dx'^) = dx^ (X'^) Aac dx' = [A^'^ Kd dx' = dx^; p^s(df) = 0, 

for the functions f G C°°(M) such that X^if) = 0. The dual projector pis: TW|m ^ TS 
is then given by 

(16) PTs =-X'^Abadx^ 
as wg is skew- symmetric. Then it follows that 

PTs(X^) =-X'^AbadxMX^) =X^; pTs(Xf) = 0, 

and so the composition {ts = tlPis- T*S TW\s gives jjsldx'^) = X'^. The Poisson 
structure TTg on W is obtained from ws, 

Hs = ^ Aab Udx''] A lls(dx^) = 1 A^b X'^ A X^ 

or directly from FTw by projection: 

ns(a,(3] =pTs(nw)(cx,|3) =nw(p^s('^),P^s((3]). 
The complementary relation gives the transverse structure TTm: 

TTm = TTw — Hs = TTw — j '^^b A X^. 
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As 

ns (df , dg) = 1 Aab X'^ A X'' (df , dg) = {f , x^jw A.b {x^ g}w, 
for extensions f, g, if we set, 

{f, g}M = nw(df, dg) - ns(df, dg), 
then it follows the well known Dirac's bracket formula: 

(17) {f , giM = {f , glw - {f , ^''Iw Aab {x^, gV- 

Note that, if we take into account (IHj) and (fT6ll . then 

ws(pTs(Xf), PTs(Xg)) = dx^(Xf) \a dx^ (Xg) A^b wg (X^ X^) 

= {f, x'^lwAcaix'', glw, 
as (obviously) it would be expected. 

4. Examples 

4.1. The Chaplygin-Caratheodory's sleigh. We first recall that the Euler-Lagrange 
equation for non-holonomic mechanical systems, using the projector method, is imple- 
mented as follows. 

Consider a mechanical system described by a Lagrangian L defined on the bundle TW 
of an n-dimensional configuration manifold W, L(z^,v^,t) — T(z^,v^,t) — V(z^,v'), where 
the z^'s are the coordinates of W and v^'s are the velocities, submitted to non-holonomic 
constraints of the form 

(18) Af(z^)v^ + B^(t) = 0, 1 < a <m. 

The stationarity condition of the action takes the well-known form 



(19) 



't2 




^9L\ 


9L 


- tl 


dt 


V9vV 


~ azi 



Sz^dt 



rt2 



Ei Sz^dt = 0, 



tl 



From (fT8|) . it follows that the virtual displacements 5z^ are not all independent, since 
they must satisfies m equations 

(20) Ai^(z^) 5z^ = 0, 

such that all the 

d /9L\ _ 9L 
^ ~ dt [dv^J ~ 9? 

cannot be set zero. 

In order to set up the projector method it is more convenient to express all the equations 
in the matricial form. Thus (fT9|) and (1201) are respectivelly written as 



ft2 

(21) 



6z dt = 

tl 



(22) A5z = 

where E, 6z are n x 1 matrices, A is the m x n matrix (A^) and the superscript t denotes 
transpose. 
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Next we shall assume that the Riemannian metric g is given by the kinetic energy 
T = (1/2]((v^)^ ■ • • + (v"^)^). We set v* = g*(»,v) for the associated co-vector. The 
constraint equation at the right of (!22|) splits the TW (with respect to g) such that 
TW = T> ® T)^, where T> is the space of all virtual displacements compatible with the 
constraints. If 5r\ is an arbitrary (linearly independent) virtual displacement, then the 
virtual displacement compatible with the constraints can be written as 

(23) 6z = P'5ti 

where P is the projector P: TW V. Substituition of ( 12^ in ( ETj) gives 

P^ 6ri dt = 

' t, 

Now as 5ri are linearly independent, we have 

(24) P E = 0, 

that is, the equation of motion compatible with the constraints must be such that, at each 
point, the Euler-Lagrange vector field is a vertical vector, while the virtual displacement 
6q is a horizontal vector, according to the D'Alembert Principle: E* 6z = 0. 
The equation flMl) can be re-written as 



(25) E = QE, 

where Q: TW T>^ is the complementary projector, defined by 

(26) Q = (A*)^G-^A, 

where G = A (A*)* is a non-sigular matrix for all mechanical systems known. Notice 



that the right-hand side of fl2^ gives the constraint forces. 

To illustrate these features we apply the projector method to obtain the equation of 
motion for the Chaplygin-Caratheodory's sleigh, a mechanical example of the Lagrange's 
Problem (we remit to Chaplygin, ref. [H], Caratheodore [10], Neimark & Fufaev, ref. [3lj 
and the paper of Koiller, ref. [22], in particular § 4.1, for further details; for a general 
geometric setting on the subject see ref. [7]). 

The sleigh consists of a rigid body supported on a horizontal plane, by three points, two 
of which slide freely (without friction) and the third which is a knife edge (or the edge of a 
cutting wheel) rigidly fixed on the body. We will consider here the special case studied by 
Caratheodory [IHlEni, in which the center of mass (cm) lies on the straight line, £, passing 
through the point of support, p, of the knife edge. The position of the body, in a fixed 
coordinate system in the horizontal plane, is determined by the coordinates (x, ij) G IR^ 
of p, and the angle 9 between the line I and the x-axis (and so the configuration manifold 
is W = IR^ X S^). So, when the body shdes, the velocity of p = p(t) can be decomposed 
in a component along the line £ 

(27) u = xcos9 -I- ij sinG 
and a component perpendicular to it 

(28) V = -r 6 tj cose - X sine 

where r is the distance from the cm to the point p and x, tj , e are the velocity coordinates 
(the dot indicate differentiation with respect to t) . The non-holonomic constraint is given 
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by the condition that the point p can move freely on the plane only in the direction along 
the line t, but not in the direction perpendicular to it, that is 

(29) v = -rQ + ycosd-xsinQ = 0. 

This condition is expressed, by the following restriction on the virtual displacements: 

(30) -r5e + cose5y-sine6x = 0. 
In this case, we have 

(31) A = (-sine, cose, -r) . 

To symplify, we set the mass of the system equal to unity and also we shall restrict the 
example to the simple case where no external forces or torques exists: so 

(32) L = T = l((x)2 + (i))^ + je^), 

where J is the momentum of the inertia about the symmetry axis through the cm. There- 
fore, we obtain g in matrix notation, 

/I 0^ 

(33) 10 

\0 J 



Using dsn]), dSU) and HMD we obtain 

sin^ e —cos e sin e 

J+r2 



J , , 

Q = — — J I —cos e sin e cos 6 

J-^rsine -J-^rcose J"^ 




Remark 2. The action of the constraint form defined by (1301) . 

a = —sin 9 dx + cos dij — r d9 

on the vector field 

X = — sm 0- — h cos 9- -— — 

ox oy J 99 

gives 

-i+r^ dcf „ 
(x(X) = — j— = G, 

and so Q = g* X (g) a. We notice also that ||X|p = g*(X, X) = G. □ 

To conclude this Lagrangian description, and to recover the equations found by Chap- 
lygin and Caratheodory ([lO], [Ml [22]) we first obtain the projected equations of motion 
taking into account ( 125|) : 

(x\ J / ^^^^ ^ —cos 9 sin 9 r sin 9 
y = J -cos 9 sin 9 cos^ 9 -r cos 9 
9/ -^^^ \J-irsin9 -J-^rcos9 J-^ 

So, taking into account the constraint equation fl30|) . its first derivative, and the system 
( IMj) . we obtain the Chaplygin-Caratheodory equations 




(35) d) = -\^j^ ) uo), u 

where cu = 9. 
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The Hamiltonization is straightforwardly obtained, since the fiber derivative of L 

is a diffeomorphism. Locally, Ci[x,y,Q,x,y^Q) — (x, y, 9, px, Py, pe), pi = (9L/9v^), 
where the v's are the x,ij,9, and the p's are (9L/9x), etc. Therefore, the constraint 
submanifold is defined by — sin9pi + cos9p2 — (r/J]pe = and the Hamiltonian is 
H = (l/2)(p^ + Py + (VJ)Pe)- Clearly, we reproduce the same arguments as above to 
obtain the constrained equations. 

4.2. A particle subjected to the constraint w = dz — ydx. Let us consider a free 
particle in M = IR^, with coordinates (x, ijjZ), subjected to the non-holonomic contact 
form w = dz — ydx. The motion of the particle in the configuration space is given by 
the Lagrangian L = [^ /2)[[v^)^ + (v^)^ + (v^)^), where {v^,v^,v^) are the corresponding 
coordinates on the fibers of TM, and the Riemannian metric g on M is given by L. So g 
is the identity tensor dx (g) dx + dy dy + dz dz, and the metric expressed in the basis 
{dx, dy, w} is 

g* = ( 1 + y ^] dx ® dx + y dx ® w + dy ® dy + y w (g) dx + w ® w. 
Notice that (see [7]) 

gaa - gap 7^ 0. 

With respect to the dual basis {Xi — 9/9x + y 9/9z, X2 = 9/9y,Y3 = 9/9z} one has, in 
matrix notation, 






1+y^ 

and G = (1 + y^) 9/9z ® 9/9z implies G(w,w] = (1 +y^) ^0. Now, 

9 9 9 
g.(w, .] = -y Xi + (1 + y^) - = - - y - = ^, 



9z 9z 



9x 



and so (IE]) is given by 

Q = g^(.,£,)[G(w,w) 



1 



l+y2 
1 

TTu2 



W ' 



[y^ dx < 



9 

9^ 



y dx' 



9 

9^ 



y dz< 



9x 



+ dz< 



4' 



Therefore, in matrix notation, Q and P are 

Q ' 



1 












e 







, p = 












1 +y^ 





1 





1+y2 


T + y^/ 



The P-vectors are 



X 



9 9 p _ 9 

and the Q- vector X*^ is £,. Again, remark that ||X^|| = 
we may also write Q — (||X^||)^^ g*(*) X^) g*(w, •). 



'(XQ,XQ) = (1 +y^) and so 
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The relationship of the above Q with the orthogonal case, that is, gaa — g«(3 = is 
the following. Let us denote by q the projector Q for this new situation (see f lTUj) ). Then 

/ 0^ 
q= 

\-v 1 

As the P- vector fields are in the kernel of q, we deduce that — ^i—^i, and 
X^i = 9/9z. 

It is easy to verify that 

[X^,Xf]=X^; [XP,X'i]=0; [X^,,X^]=0. 

This is just the Heisenberg algebra, a very well known non-holonomic Lie algebra (and a 
fundamental example in sub-Riemannian geometry). Also, 

/I +1)2 -y' 
U= 10 
\ 1 

is so that Q = U q U"^ . 

Let us now consider the following constraint equation 



(36) f(v%vy,v^ 
for this problem. The use of 



ijv" = 0. 



(37) 



gives 



l+ij2 






1 





1+r 





yz]=0, i) = 0. 



Therefore 
(38) 



x + -yz = 0, y = 0. 

Equation 0381) is just the momentum equation of Bloch et Al. |3] (page 85). We remark 
that from (15^ one obtains z — yx + yx and so, using this equation in ( 155]) we obtain the 
Bates-Sniatycki motion's equations (ref. |1]) 

V 



x + 



1 



xy 



0, 



0. 



The motion of the particle in the phase space (endowed with the natural Poisson struc- 
ture TT) is given by the Hamiltonian 

restricted to the submanifold defined by the equation — yp^ = 0. The momenta, 
compatible with the constraint, can be written as p = P p or = P|pj, where i, j = 
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x,y,z. Thus, 
(39) 




/_L_ 








\1 



1 

,.2 

T 




Note that = (px + yPz) = yPx- The apphcation of the local definition of the 
Poisson brackets on a regular Poisson manifold gives 



{z\z^} = 0; K,pJ = P|; {p„pj 




and so, after using the constraint equation p^ = ypx, we obtain the following pseudo- 
Poisson structure, that is, the Poisson bracket is skew-symmetric, satisfies de Leibniz rule 
but it may not satisfy de Jacobi identity (see Marie [30], for further details): 



(40) 



n 








1 





















V 1 







1 +-y^ 


,,2 



1 



1 










1 









Finally, we want to remark that there are many others non-holonomic systems which 
admits a pseudo - Poisson structure. The Chaplygin-Caratheodory's sleigh (example 4.1) 
is the protoptype of mechanical systems described by the Poisson geometry. 

In fact, most of the known non-holonomic mechanical systems are given by completely 
non-integrable constraints which are linear in the velocities (as in the examples 4.1 and 
4.2). In these cases, as we have shown in 4.2, the induced Poisson structure is given by: 



(41) 



n 



P 

-p D 



where P is the n x n matrix-projector and D is the n x n matrix whose elements are 

9P^ 





Generally, the matrix TT is singular because the matrix P is generally singular, that is 
det P = . We say "generally" because one could conceive cases in which the system 
of constraints is completely non-integrable but - nevertheless - admits a singular integral 
which reduces the system to the integrable form. In this case, as in the holonomic case, 
TT, is invertible, which gives rise a symplectic structure. The authors do not know any 
concrete examples of non-holonomic mechanical systems of this kind. 
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